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The interaction of two rectangular pulses in nonlinear Schrodinger model is studied by solving the 
appropriate Zakharov-Shabat system. It is shown that two real pulses may result in appearance of 
moving solitons. Different limiting cases, such as a single pulse with a phase jump, a single chirped 
pulse, in-phase and out-of-phase pulses, and pulses with frequency separation, are analyzed. The 

ff) . thresholds of creation of new solitons and multi-soliton states are found. 
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I. INTRODUCTION 



The nonlinear Schrodinger (NLS) equation is an important model of the theory of modulational waves. It describes 
i— i . the propagation of pulses in optical fibers 1,2 , the dynamics of laser beams in a Kerr media, or the nonlinear difraction 3 , 
waves in plasma 4 , and the evolution of a Bose-Einstein condensate wave function 5 . The NLS equation is written in 
dimensionless form as the following: 

iu z + u xx /2 + \u\ 2 u = 0, (1) 

where u(x, z) is a slowly varying wave envelope, z is the evolutional variable, and x is associated with the spatial 
variable. 

An exact solution of the NLS equation has a form of a soliton: 

u(x, z) — 2i] sech[2^(a; + 2(z — xq)] exp[— 2i£x — 2i(£ 2 — -q 2 )z + i^o] • (2) 

O ■ 

where 2r\ and 2£ are amplitude, or the inverse width, and the velocity of the soliton, xq and 4>o are the initial position 
and phase, respectively. The soliton represents a basic mode and plays a fundamental role in nonlinear processes. 
The dynamics of NLS solitons and single pulses even in the presence of various perturbations is well understood 
(see e.g. Refs. 1,2 and reference therein). However, the evolution of several pulses is not studied in details. In 
recent works 6 (see also book 2 ) mostly an interaction of solitons and near-soliton pulses was considered. A study 
of near-soliton pulses, especially a use of the effective particle approach, often results in small variation of soliton 
parameters, including the soliton velocities and as a consequence weak repulsion or attraction of solitons. Such a 
study do not involve a possibility of appearance of additional solitons. However, for many applications it is necessary 
to consider the interaction of pulses with arbitrary amplitudes or pulses with different parameters. For example, 
in optical communication systems with the wavelength division multiplexing (WDM), the initial signal consist of 



several solitons with different frequencies. An estimation of the critical separation between pulses is important for 
- - 1 determination of the repetition rate of a particular transmission scheme. 

In present work the interaction of two pulses in the NLS model is studied both theoretically and numerically. We 
show the presence of different scenarios of the behaviour, depending on the initial parameters of the pulses, such as 
the pulse areas, the relative phase shift, the spatial and frequency separations. One of our main observation is a fact 
that a pure real initial condition of the NLS equation can result in additional moving solitons. As a consequence the 
number of solitons, emerging from two pulses separated by some distance, can be larger than the sum of the numbers 
of solitons, emerging from each pulse. Such properties were also found for the Manakov system 7 , which is a vector 
generalization of the NLS equation. The scalar NLS equation was studied in 7 as a particular case. Later similar 
results and approximation formulas for the soliton parameters were obtained in papers 8 ' 9 (see also 10 ). In works 
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mostly the interaction of real pulses was analyzed, while here we consider pulses with non-zero relative phase shift 
and frequency separation. A preliminary version of this study was presented in work 11 . 

The paper is organized as the following. The linear scattering problem associated with the NLS equation is 
considered in Section II. We also present the general solution of the problem for the case of two rectangular pulses. 
In Section III we study different particular cases, such as two in-phase pulses, two out-of-phase pulses, a single pulse 
with a phase jump, a single chirped pulse, and two pulses with the frequency separation. The results and conclusions 
are summarized in Section IV. 



II. DIRECT SCATTERING PROBLEM 



In this paper we are interested only in an asymptotic state of the pulse interaction. In order to simplify the problem 
and to obtain exact results we consider the interaction of two rectangular pulses ("boxes"). Therefore we take the 
following initial conditions for Eq. (1): 



Qi exp[2iuix], for x\ < x < X2, 
u(x,0) = U(x) — ^ Q 2 exp[2iu 2 x], for x 3 < x < x&, 
0, otherwise , 



(3) 



where Qi and Q 2 are the complex constant amplitudes, w\ = x 2 — x\ and w 2 = X4 — x 3 are the pulse widths, and 2v\ 
and 2v 2 are the detunings. 

It is known that NLS equation is integrable by the inverse scattering transform method 3 . As follows from this fact, 
initial conditions, which decreases sufficiently fast at x = ±00, results in a set of solitons and linear waves (so called, 
radiation). The number N and parameters of solitons emerging from an initial condition are found from the solution 
of the Zakharov-Shabat scattering problem 3 : 



-i^-iU*{x)^=\i, 2 



with the following boundary conditions 



i\x 



/a(A), 
\b(X)e 



Ax 



(4) 



(5) 



Here ^(x) is an eigenvector, A is an eigenvalue, a(A) and 6(A) are the scattering coefficients, and a star means a 
complex conjugate. The number N is equal to the number of poles A„ = £„ + ijj n , where n — 1, . . . , N, and r\ n > 0, of 
the transmission coefficient l/a(A). Each A„ is invariant on z. If all are different then u(x, z) at z — > 00 represents 
a set of solitons, each in the form of Eq. (2) with r/ = t] n and £ = £„. If real parts of several A„ are equal then a 
formation of a neutrally stable bound state of solitons is possible. 

The solution of the Zakharov-Shabat problem (4) with the potential (3) is written as 



a(A) = e 



n \ -( A + ^l) • (i \ 

cos(fc 1 ui 1 ) — 1 sm(fciwi) 

ki 

11 \ -(A + ^2) . . 
cos{k 2 w 2 ) — i ; sm(fc 2 W2) 



k 2 

QXQl c _2»(A+^i)x2 e 2i(X+u 2 )x 3 

k\k 2 



sin(fciwi) sin(fc 2 W2) 



(6) 



6(A) 



i(A+L/i)«;i — i{\+v 2 )w 2 



(_Ql e -2i(X+, 1 ) X2 sin(fclU;i) 
I fcl 

n n . .(A + ^2) . n , 

cos(k 2 w 2 ) + i sm(fc 2 W2) 

k 2 



2 



Q2 e -2i(X+, 2 ) X3 sin(fc2W2 ) 
^2 



COS 



:(kiWi) 



■ (A + 1/1) 



fcl 



sin(fciwi) 



(7) 



where ^ = [(A + ^) 2 + |Q,f] 1/2 - 

Since the linear operator in (4) is not Hermitian, complex eigenvalues are possible even for real u(x, 0) (e. g. 
see Section III A). Though this is an obvious fact, "an interesting "folklore" property seems to have arisen in the 
literature over the last 25 years, namely, that only pure imaginary EVs (eigenvalues) can occur for symmetric real 
valued potentials" 8 . As demonstrated below, the statement [Theorem (III) in §2] of paper 12 , that claims this result, is 
incorrect. An existence of eigenvalues with non-zero real parts for Zakharov-Shabat problem with pure real potential 
was first shown in paper 7 . 

Equations (6,7) represent a general solution of the scattering problem (4,5) with initial condition (3). Applications 
of these equations to particular cases of the pulse interaction are considered in the next section. 



III. RESULTS 



A. Interaction of in-phase pulses with equal amplitudes 

1. Properties of eigenvalues 

Here we analyze a simple case of two real pulses, separated by a distance L = x 3 — x 2 , with zero detuning, i.e. 
Qi = Qi = Qo, w i = w 2 = w, and v\ = u 2 = 0, where Qo is real. Then using Eq. (6), the equation for discrete 
spectrum is written as 

F(X,Q ,w)±^-e iXL sm(kw) =0 , (8) 
k 

where F(X, Q, w) = cos(fc w) — iAsin(fc w)/k, and k = (A 2 +Q 2 ) 1 / 2 . Note that F(X, Qo,w) = determines the discrete 
spectrum for a single box with zero detuning 13 . Therefore the second term in Eq. (8) can be associated with the 
result of nonlinear interference. Recall also that for a single box with amplitude Qo and width w, the number N$b 
of emerging solitons is determined as 3 N$b = int(Qow/ir + 1/2), where intQ means an integer part. Results for the 
two boxes are reduced to those for a single box in limiting cases L = and L = oo. 

As shown by Klaus and Shaw 8 , the Zakharov-Shabat problem with a "single-hump" real initial condition admits 
pure imaginary eigenvalues only, i.e. solitons with zero velocity. We show that the case of two pulses provides much 
richer dynamics. 

Let us now compare the properties of eigenvalues at different S = Qqw (Fig. 1). In Fig. 1, as well as in subsequent 
figures of the paper, all variables are dimensionless. In the first two cases, S = 1.8 and S = 2.0, there is one soliton 
at L — and there are two solitons at L — oo, while in the case S — 2.5 there are two solitons in both limits. The 
dependence of eigenvalues on L at S = 2.5 is obvious, while that at S = 1.8 and 2.0 looks unexpected. Firstly, the 
number of solitons at intermediate L is larger then that in the limits L = and L = oo. Secondly, the two real boxes 
lead to eigenvalues with non-zero real part. Third, for S = 2.0 there is a "fork" bifurcation at L — Lp w 4.1, when 
two eigenvalues coincide. At larger L three pure imaginary eigenvalues constitute a three-soliton state, so that the 
limiting two-soliton case at L — > oo is realized as a limit of a three-soliton solution with an amplitude of the third 
soliton tending to zero. 

Results of numerical simulations of the NLS equation (1) agree with analysis of Eq. (8). For example, as shown 
in Fig. 2, in accordance with Fig. lb there are one fixed and two moving solitons at S = 2.0 and L = 2, and there 
are a three-soliton state and two moving solitons at S = 2.0 and L — 5. Note that an appearance of moving solitons 
and multi-soliton states is not related to the rectangular form of initial pulses. For example, an initial condition 
u(x, 0) = 0.7[sech(x + 2.5) + sech(x — 2.5)] also results in moving solitons. 

Below we discuss in details the behaviour of the eigenvalues, namely we find a threshold of appearance of new roots, 
estimate a number of emerging solitons, and calculate a threshold for the "fork" bifurcation. It should be mentioned 
that eigenvalues with non-zero real part do not exist only at S = [37r/4, 3.3] and S = [7ir/4, 5.51] (see Section III A2), 
so that the dependence at S = 2.5 is rather an exception than a general rule. This results allows to understand why 
moving solitons are not observed in interaction of near-soliton pulses with area S w w. 
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2. Appearance of new eigenvalues 



Solving numerically Eq. (8), one can conclude that new eigenvalues penetrate to the upper half-plane of A in pairs 
by crossing the real axis. Therefore, the bifurcation parameter can be found from Eq. (8), assuming that A = (3, 
where (3 is real: 

\/2S 2 — v 2 

cotan y = ±-¥ — — , (9) 

(3 = ±Q sm(f3L) . (10) 

Here y = kw, k — (/3 2 + <3o) 1/2 , and the signs are taken such that tan(y) tan(/3L) < is satisfied. As follows from 
the definition of y and Eq. (9), one has S < y < 2S. 

Analysis of Eqs. (9,10) results in the following conclusions: 

(i) As follows from Eq. (9), the number Npp of the penetration points depends only on S and is determined from: 

Np P = 4(m-ra+l) -26 
3 
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n + - 



A6(S m - S) for S > 3tt/4 , (11) 



where to = int(v / 2 S/ir), n = int(5/7r), 6(x) is the Heavysidc function, and S m is a root of 

tan(v/25^ - 1) = ^2S 2 n - 1 , (12) 

which satisfies mix < (25^ - l) 1 / 2 < (to + l)ir. It is easy to find that N PP = for S < ir/4 and N PP = 2 for 
7r/4 < S < 37r/4. Equation (12) defines such values of S = S m , when the right hand side of Eq. (9) with plus sign 
touches cotany curve. All penetration points (3j, where j — 1, . . . ,Npp, are symmetrically situated with respect to 
[3 = 0. 

(ii) All roots \(3j\ < Q , which follows from 2S 2 - y 2 > 0. 

(iii) For every (3j, Eq. (10) defines the separation distance L = Lq, when eigenvalues cross the real axis. 

(iv) As follows from Eq. (10) there is an infinite number of thresholds Lq for a given (3j. However, the total number 
of eigenvalues in the upper half-plane of A is, most probably, finite, because for some Lc eigenvalues pass to the upper 
half-plane, and for other Lc eigenvalues go to the lower half-plane. The direction of eigenvalue motion is defined by 
the derivative dX/dL at A = (3j. 

The positions of penetration points, (3j as a function of S is shown in Fig. 3a, where only positive j3j arc presented. 
As follows from Eq. (9) the number Npp decreases by two, when S passes (21 + l)7r/4, where I = 1,2..., and Npp 
increases by four, when S exceeds S m [see Eq. (12)]. Therefore one can obtain that Eq. (9) has no roots only at 
S = [37r/4, S2) and at S = [7n/4, S3], where S2 ~ 3.26 and S3 w 5.51 are found from Eq. (12). This property is clearly 
seen in Fig. 3. The dependence of Lc on S is presented in Fig. 3b. Only the thresholds, such that (3jLc = [0, 2w], 
are shown for each (3j. 



3. Thresholds of "fork" bifurcation 

Here we analyze a bifurcation, when a pair of complex eigenvalues becomes pure imaginary, e.g. Lp w 4.1 in Fig. lb. 
The equation which determine pure imaginary eigenvalues can be obtained from Eq. (6) with Re [A] = 0, i.e. A = 27: 

cotan „ = -VS^^SeM-VS^LM 

V 

Here y — kw , k = (— 7 2 + Qq) 1 ^ 2 . It is easy to show that k 2 should be positive (there is no real solution for k 2 < 0). 
As a consequence, all pure imaginary eigenvalues satisfy 7j < Qq. 

The value of L = Lp, when new pure imaginary root of Eq. (13) appears, corresponds to the "fork" bifurcation. 
The bifurcation threshold Lp can be found from the condition that the functions, corresponding to the right-hand 
side of Eq. (13), touch the cotan y curve. 

Low bound of the number of the pure imaginary eigenvalues is found as Nlb = int(2S'/7r + 1/2). This number is 
an actual number of the pure imaginary eigenvalues for all S, except of the regions where 
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- + irl < S < — + 1 = 0,1,. 
2 4 



(14) 



If S satisfies Eq. (14) then the number of pure imaginary eigenvalues can be either Nlb or Nlb + 2, depending on 
whether L < Lf(S) or L > Lf(S), respectively. Therefore the appearance of new pure imaginary eigenvalues, in 
other words, the fork bifurcation, is possible only if 5* satisfies Eq. (14). Figure 4 represents the dependence of hp on 
S, where only one interval, corresponding to I = 1 in Eq. (14), is shown; the behaviour for I > 1 is similar. One can 
calculate that L F (S = 1.8) = 10.4, that is why the fork bifurcation is not seen in Fig. la. 

B. Two out-of-phase pulses with equal amplitudes 

In this section we study the influence of constant phase shift on the pulse interaction, i.e. we consider Qi = 
Qoexp(—i a), Q2 = Qoexp(i a), where Qo and a are real, Wi = W2 = w, and v\ = V2 = 0. The non-zero relative 
phase shift, 2a, changes greatly the properties of the eigenvalues, so that the behaviour presented in Section III A is 
hard to realize in experiments, because it is difficult to prepare two pulses exactly in phase. The phase shift breaks 
the simultaneous appearance of a pair of solitons at A = ±/3j, and affects to the "fork" bifurcation. 

For a^O, the equations for eigenvalues and for penetration points can be obtained from Eq. (8) and Eqs. (9,10) 
by changing XL — ► XL + a and (3L — ► (3L + a in the exponent and sinus functions, respectively. Therefore the number 
and positions of the penetration points are the same as for the case a = 0. As for the threshold Lc, it is shifted on 
the value a/(3j, so that Lc{a) = Lc(a = 0) + a/0j, where only Lc > should be taken into account. 

The influence of the phase shift on the distribution of eigenvalues is shown in Fig. 5. As seen, now new eigenvalues 
appear one by one, not in pairs, and, as a concequense, the "fork" bifurcation disappears. Further, the real parts of 
the roots do not vanish at finite L, but decreases smoothly. This means that the presence of the phase shift breaks 
up a multi-soliton state, which is known to be neutrally stable to perturbations. 

At L = 0, the phase shift corresponds to the phase jump of a single pulse. Such a phase jump can result in an 
appearance of additional solitons as shown in Fig. 5b. The threshold of the phase shift a t h, when the first new soliton 
appear, can be found from the condition a t h = \Pi\Lc{a = 0), where (3\ is the position of the penetration point 
nearest to zero. 



C. Two pulses with frequency separation 

In this section we analyze the initial condition (3) with the following parameters Qi = Q2 = Qo, »i = W2 = w , 
—v\ = V2 = v, where Qo is a real constant. This case models the wavelength division multiplexing in optical fibers, 
the case when an input signal consists of two or more pulses with different frequencies. Actually, since Qo can be taken 
sufficiently large we consider the interaction of multi-soliton states. The detailed analysis of the interaction of sech- 
pulses at different frequencies is presented in papers 14 and in review 15 . In particular, the authors of papers 14 ' 15 consider 
the evolution of a superposition of N solitons with the same position of the centers, but with different frequencies. As 
shown in these works there is a critical frequency separation, above which N solitons with almost equal amplitudes 
emerge. Below this critical value the number of emerging solitons can be not equal N and their amplitudes can 
appreciable differ from each other. It was also demonstrated that an introduction of a time shift between pulses 
results in decrease of the frequency separation threshold. The geometry described by Eq. (3) corresponds to the 
combination of WDM and time-division multiplexing schemes, therefore our study can give some insight to such a 
behaviour of the threshold. Moreover, the authors of works 14,15 mostly used the perturbation technique and numerical 
simulations, while in the present paper we deal with an exact solution of the Zakharov-Shabat problem. 

First let us consider the case L = that correspond to the case of a single chirped pulse of width 2w. The 
dependence of the eigenvalues on v, which plays here the role of a chirp parameter, is shown in Fig. 6. At small v the 
interaction of the pulse components is strong, so that there is one pure imaginary eigenvalue, or a single soliton with 
zero velocity. At larger v the frequency difference of the pulse components results in a repulsion of the components, 
or a pulse splitting. At sufficiently large v the velocities of emerging solitons tends, as expected, to ±2za 

There is also a narrow region of v, e.g. v = [0.98,0.99] on Fig. 6, where three solitons, one fixed and two moving 
solitons, exist. This region separate two different types of the evolution of a chirped pulse. The left boundary, which 
corresponds to the appearance of new eigenvalues, of the region is found from the condition similar to that considered 
in Section III A 2. The right boundary is found from a(X = 0) = 0, which defines the values of v as a function of the 
other parameters, when the pure imaginary root disappears. 

The dependence of the spectrum on L is presented in Fig. 7. At small v (Fig. 7a) we see again an appearance of 
additional solitons similar to the case v = (Fig. 1). At larger v (Fig. 7b) the repulsion is so strong that it suppresses 
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the appearance of small-amplitude solitons. Therefore there is a threshold of the frequency separation above which 
the interaction of two pulses is negligible. This result is in agreement with the conclusions of the paper 15 . 



IV. CONCLUSION 



The interaction of two pulses in the NLS model is studied by means of the solution of the associated scattering 
problem. The strong dependence of the dynamics on the parameters of the initial pulses is shown. For intermediate 
separation distances L the existence of additional moving solitons is possible even in the case of two in-phase pulses with 
the same frequencies. These additional solitons can be considered as a result of nonlinear interference of pulses. The 
phase shift of two pulses removes a degeneracy in the behaviour, namely it affects to the symmetry of the parameters 
of emerging solitons and results in a break-up of multi-soliton states peculiar to the in-phase case. It is also shown 
that the strong frequency separation suppress the appearance of additional solitons. The results obtained in the 
present paper can be useful for analysis of the transmission capacity of communication systems and for interpretation 
of experiments on the interaction of two laser beams in nonlinear media. Recently, the generation of up to ten solitons 
has been observed experimentally in quasi-lD Bose-Einstcin condensate of 7 Li with attractive interaction 16 . Our 
study can be also helpful for interpretation of this experiment. 



ACKNOWLEDGMENTS 



This researh was partially supported by the Foundation for Support of Fundamental Studies, Uzbekistan (grant N 
15-02) and by FAPESP (Brasil). 



1 G. P. Agrawal, Nonlinear Optics (Academic Press, San Diego, 1989); F. Kh. Abdullaev, S. A. Darmanyan, and P. K. 
Khabibullaev, Optical Solitons (Springer- Verlag, Heidelberg, 1993). 

2 A. Hasegawa and Yu. Kodama, Solitons in Optical Communications (Clarendon Press, Oxford, 1995). 

3 V. E. Zakharov and A. B. Shabat, Sov. Phys. JETP 34, 62 (1972). 

4 See e. g. H. Ikezi, In Solitons in action, edited by K. Longren and A. Scott, (Academic Press, New York, 1978). 

5 E. P. Gross, Nuovo Cimento 20, 454 (1961); J. Math. Phys. 4, 195 (1963) L. P. Pitaevskii, Zh. Eksp. Teor. Fiz. 40, 646 
(1961) [Sov. Phys. JETP 13, 451 (1961)]. 

6 C Desem and P. L. Chu, In Optical Solitons - Theory and Experiment, edited by J. R. Taylor (Cambridge Univ. Press, 1992), 
Chap. 5; D. Anderson and M. Lisak, Opt. Lett. 18 , 790 (1986); V. I. Karpman and V. V. Solov'ev, Physica D 3, 487 (1981). 

7 F. Kh. Abdullaev and E. N. Tsoy, Physica D 161, 67 (2002). 

8 M. Klaus and J. K. Shaw, Phys.Rev E. 65, 036607 (2002). 

9 M. Desaix, D. Anderson, L. Helczynski, M. Lisak, In Nonlinear Guided Waves and Their Applications, technical dijest of 
the OSA International Workshop (Streza, Italy, 2002), vol.80, NLTuD13. 

10 After the submission of this paper we became aware with the work by M. Desaix, D. Anderson, L. Helczynski, and M. 
Lisak, Phys. Rev. Lett., 90, 013901 (2003), which deals with real initial conditions. However, here we consider more general 
complex initial conditions (see Eq. (3)) and study different scenarios of the pulse interaction depending on space and frequency 
separations, phase shift and pulse areas. 

11 E. N. Tsoy and F. Kh. Abdullaev, In Nonlinear Guided Waves and Their Applications, technical dijest of the OSA Interna- 
tional Workshop (Streza, Italy, 2002), vol.80, NLTuD14. 

12 J. Satsuma and N Yajima, Suppl. Prog. Theor. Phys. 55, 284 (1974). 

13 S. V. Manakov, Sov. Phys. JETP 38, 693 (1974). 

14 P. A. Andrekson, N. A. Olson, J. R. Simpson, T. Tanbun-Ek, R. A.Logan, P. C. Becker, and K. W. Wecht, Appl. Phys. Lett. 
57, 1715 (1990); Y. Kodama and A. Hasegawa, Opt. Lett. 16, 208 (1991); C. Etrich, N.- C. Panoiu, D. Mihalache, and F. 
Lederer, Phys. Rev. E. 63, 016609 (2000). 

15 N. -C. Panoiu, I. V. Mel'nikov, D. Mihalache, C. Etrich, and F. Lederer, J. Opt. B: Quantum Semiclass. Opt. 4, R53 (2002). 

16 K. Strecker, G. Partridge, A. Truscott, and R. Hulet, Nature 417, 150, (2002). 



G 



012345678 
L 



01 2345678 
L 




' 1 1 1 1 1 1 1 1 

012345678 

L 

FIG. 1. In-phase pulses: The dependence of real (dashed lines) and imaginary (solid lines) parts of A n on the separation 
distance, w = 1. Numbers near lines corresponds to n. (a) Q — 1.8, (b) Q = 2.0, (c) Q = 2.5. 




FIG. 2. Evolution of two rectangular pulses, Qo = 2,w = 1. (a) One fixed soliton and two moving solitons at L — 2. (b) 
Three-soliton state at L = 5. 
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FIG. 3. (a) The dependence of f3j on S. (b) Threshold Lc, when eigenvalues cross the real axis of A-plane, as a function of 
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FIG. 4. Threshold Lf of the fork bifurcation as a function of S. 
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FIG. 5. Out-of-phase pulses: The dependence of real (dashed lines) and imaginary (solid lines) parts of \„ on L for 
Qo = 2, w = 1. Numbers near lines corresponds to n. (a) a — n/8, (b) a = 7r/4. 
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FIG. 6. Single chirped pulse; The dependence of real (dashed lines) and imaginary (solid lines) of A n on v for 
Qo = 2, w — 1, L — 0. Numbers near lines corresponds to n. 
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FIG. 7. Pulses with frequency separation: The dependence of real (dashed lines) and imaginary (solid lines) of \„ on L for 
Qo = 2,w = 1. Numbers near lines corresponds to n. (a) v — 0.5; (b) v = 1.25. 
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